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Abstract 


Rb3;D(SeO,), belongs to the well-known group of hydrogen sulfate and selenate crystals which are promising for fuel cell applications. However, 
the high temperature properties of this salt have been much less extensively studied than those of Rb3H(SeO,)2. Superionic phase transition in 
Rb3D(SeO,). and Rb3H(SeO,)2 single crystals was studied using impedance spectroscopy and DSC methods. Temperature evolution of the 
ferroelastic domain structure was observed under a polarizing microscope. Additionally, the X-ray diffraction and Raman scattering measurements 
for both crystals were carried out at room temperature. Rb3D(SeO,)2 undergoes a structural phase transition from the low temperature ferroelastic 
phase to the superionic, paraelastic one. A correlation was found between the ferroelastic domains evolution, the anomalies of DSC and the 
temperature dependence of ionic conductivity. The temperature range of superionic phase transition in the first heating run was much wider 
and the Ts value was lower in Rb3D(SeOx), than in the non-deuterated compound. Additionally, a considerable temperature hysteresis of Ts in 
Rb3D(SeO,)2 was observed, unlike in Rb;H(SeO,),. The activation energies of H* and D* conductivity along the c-axis were found to be almost 
the same. The temperature variation of ionic conductivity in Rb;D(SeO,), obtained in the first heating run of a virgin sample was found to differ 
considerably from that measured in subsequent cooling/heating runs. The small H/D isotope effect of proton conductivity in the crystals studied 


pointed to a classical activated diffusion mechanism. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


A common feature of hydrogen sulfate and selenate crys- 
tals of the general chemical formula A,Hy(XO4)(4)/2, where 
A denotes Rb, Cs, NH4, and X is S or Se, is the hydro- 
gen bonding. The structural motif of the crystals consists of 
XO; tetrahedra linked with short hydrogen bonds into isolated 
dimers (in the group of x=3 and y=1), into infinite chains 
(when x=1 and y=1) and into trimers (for x=4 and y=2). At 
the superionic phase transition temperature Ts, ranging from 
300K for (NH4)3H(SeO4)2 to ~490K for Rb3H(SO4)2, the 
materials undergo a transition to the phase with a fast proton 
transport. In some compounds the temperature Ts coincides 
with the ferroelastic—paraelastic phase transition temperature 
Tc. 

Although all the above crystals are water soluble and mechan- 
ically “plastic” in the superionic phase, a successful operation of 
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a fuel cell with CsHSO, (Ts = 141 °C) as a membrane material 
was reported a few years ago [1]. Acid salt electrolyte-based fuel 
cell offers several advantages over a polymeric one. The salts 
are “true” proton conductors, with no water involved in the con- 
duction mechanism, so careful water management in a fuel cell 
is not necessary. Furthermore they exhibit superprotonic con- 
ductivity at elevated temperatures, thus the purity requirements 
on the hydrogen fuel are relaxed, catalyst activity is increased 
and waste heat is more easily manageable. At present a major 
shortcoming of hydrogen sulfate and selenate electrolytes is their 
catalyzed reduction under hydrogen atmosphere, accompanied 
by the generation of H2S, which is an effective catalyst poison 
[2]. 

Rb3H(SeO4)2 and Rb3D(SeO4)2 crystals are ferroelastic, 
with C2/c space group at room temperature and transform 
to superprotonic (paraelastic) phase of R3m space group at 
Ts. In the low temperature phase the structural motif con- 
sists of isolated dimers formed by asymmetric SeOq tetrahedra 
linked with short H/D bonds between apical O(2) oxygen 
atoms of neighboring molecules. Protons (deuterons) occupy 
one of the two potential minima between their respective donor 
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and acceptor, so the bonds are ordered (i.e. their direction 
in the lattice does not change). The pattern is referred to as 
the so-called zero-dimensional hydrogen-bond network. The 
O(2) atoms and therefore the H-bonds form layers parallel 
to the (001) plane [3]. In this phase the number of struc- 
turally equivalent positions for protons equals their number, 
so any proton conductivity has to involve defect sites. As a 
result, the conductivity is low and its activation energy value 
high. 

Above Ts, in the trigonal symmetry phase, two additional 
directions of hydrogen bonds in the (001) plane are symme- 
try equivalent. The directions differ by +27r/3 and —27/3 from 
that in the monoclinic phase. In this way the number of struc- 
turally equivalent positions for protons (deuterons) increases 
three times [3]. The proton transport proceeds via structure diffu- 
sion (Grotthus mechanism) and consists of two steps: (1) proton 
transfer in an H-bond, from donor to acceptor and (2) reori- 
entation of the HSeO4 group (which involves breaking of the 
H-bond) to a position at which a new bond can be formed. It 
means that a proton encounters two energy barriers along its 
conduction path: the one separating the energy minima for trans- 
fer from donor to acceptor and the barrier for the HSeO4 group 
reorientation. 

The fast proton transport in the crystals is conditioned by 
molecular dynamics of the proton environment which involves 
both vibrational and librational modes of the SeO4 tetrahe- 
dra and lowers the two barriers for proton transport [4]. As 
a result, each proton forms H-bonds with three O(2) oxygen 
atoms of the nearest SeO4 groups with the same probabil- 
ity p=1/3. In this way protons (and consequently H-bonds) 
are de-localized in the superionic phase and form a dynam- 
ically disordered network of H-bonds in the (001) plane 
[5,6]. The arrangement of H-bonds in layers results in a dis- 
tinct anisotropy of proton conductivity, which is ~20 times 
higher in the plane than along the direction perpendicular to it 
[T7]. 

At low temperatures a dramatic H/D isotope effect in the title 
compound was reported. Rb3D(Se0O4)2 undergoes antiferroelec- 
tric phase transition at 92 K, whereas the protonated crystal does 
not exhibit any transition below room temperature [8,11,12]. On 
the other hand, insignificant isotope effect on the superionic tran- 
sition was observed in other acid salt CsHSO4: Ts in CsDSO4 
was by ~2 K lower and the activation energy of deuteron trans- 
port by 0.01—0.02 eV higher than in CsHSOg, whereas the effect 
on the ionic conductivity was expressed by o4/a? =2.3 and 
~1.5 [13,14]. Considering a much pronounced isotope effect in 
Rb3H(SeO 4). at low temperatures, it seems interesting to study 
the isotope exchange influence on the superionic phase transition 
in this crystal. 

Impedance spectroscopy conductivity measurements of 
Rb3H(SeO4)2 have shown that the temperature Ts ranges from 
449 K to 456 K [3,8,9] and the activation energy in the high con- 
ducting phase is E4 = Ep © 0.28 eV, Ec © 0.48 eV along the a, b, 
and c crystallographic axes, respectively [7]. For Rb3D(SeOq)2, 
the values of Ts = 460 K and E; = 0.28 eV in the superionic phase 
were obtained from the conductivity measurements at 100 kHz 
[10]. 


2. Experimental 


Rb3D(SeO4)2 was obtained by 6-fold recrystallization from 
D20. Single crystals of pseudo-hexagonal habit were grown 
from stoichiometric D2O solution of the deuterated salt by 
slow evaporation at room temperature. The degree of deuter- 
ation was controlled by NMR and amounted to ~90%. 
Rb3H(SeO4)2 crystals were grown by slow evaporation at 
room temperature from aqueous solution of the salt. Differ- 
ential scanning calorimetry (DSC) studies were performed 
with Netzsch DSC-200 calorimeter at a heating/cooling rate of 
2K min™!. 

The temperature evolution of the ferroelastic domain pattern 
was studied in the (00 1) plane using a polarizing optical micro- 
scope equipped with a LINKAM cooling-heating stage. The 
temperature was changed at a rate of 0.5 K min™!. 

The ionic conductivity along the c-axis, oc, of both crys- 
tals was measured in air by the impedance spectroscopy method 
with a computer controlled HP-4192A Impedance Analyzer. The 
measuring voltage was 0.5 V, in 100-10 MHz frequency range. 
For conductivity measurements samples of ~10 mm? in area 
and ~0.4mm thick were used. The electrodes were made by 
covering the opposite (00 1) faces with silver paste. Each com- 
plex impedance plot was obtained at a temperature stabilized to 
within 0.03 K. 

NIR Raman spectra (A= 1064nm) were recorded at room 
temperature with a BRUKER IFS66 FRA106 spectrometer in 
the 180° geometry. Integrated X-ray intensities were collected 
on an Xcalibur (Oxford Diffraction Company) diffractometer 
equipped with a CCD detector, using graphite monochromated 
Mo Ka radiation. Unshaped crystals of 0.23 mm (Rb3D(SeO4)2) 
and 0.26 mm (Rb3H(SeO,4)2) in average diameter were used. 
The measurements were carried out at 15K and 295K for 
Rb3D(SeQO4)2 and 120 K and 295 K for Rb3 H(SeO.)2. Low tem- 
perature was maintained with a HeliJet gas-flow cooling system 
(Oxford Diffraction Controller), AT=0.3 K. The CrysAlis soft- 
ware version 1.170.2 (Oxford Diffraction) was used for the data 
processing. The crystal structure of both salts was determined 
and refined at all temperature points using the SHELXL-97 
program package. 


3. Results 
3.1. DSC study 


DSC measurements (Fig. 1) show that the superionic 
phase transition in Rb3D(SeO4)2 is stretched over ~18 K 
and is completed at ~451K on the first heating, whereas 
on cooling it is finished at ~430K. The transition enthalpy 
AH ~ 3.4kJ mole~!. It should be added that the temperature 
of the DSC maximum on heating is close to the temperature at 
which new twin walls start to appear on the domain structure 
observations. 

DSC spectra of Rb3 H(SeOq)2 are shown in Fig. 2. The supe- 
rionic phase transition is stretched over 6-8 K and the transition 
enthalpy AH ~3.8kJ mole~!. The transition is completed at 
~460 K on heating and at ~450 K on cooling. 
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Fig. 1. DSC scans for Rb3D(SeO4)2 obtained in the first heating/cooling cycle. 
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Fig. 2. DSC scans for Rb3H(SeO4)2 obtained in the first heating/cooling cycle. 
3.2. Impedance spectroscopy study of H* (D* ) conductivity 
Proton conductivity studies have also revealed differences in 


the behavior of Rb3 H(SeO,)2 crystal and its deuterated analogue 
(Fig. 3). The o-(1/T) dependence, obtained in the first heating run 
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Fig. 3. Temperature variation of H*(D*) conductivity for Rb3H(SeO4)2 and 
Rb3D(SeO,4)2 in the first heating run. 
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Fig. 4. Temperature variation of H*(D*) conductivity for Rb3H(SeO4)2 and 
Rb3D(SeO,)2 in the first cooling run. 


for Rb3D(SeO4)2 crystal along the c-axis, differs considerably 
from the dependences measured in subsequent cooling/heating 
cycles. The distinct conductivity change at Ts, characteristic of 
the superionic phase transition [15], appears at a lower tem- 
perature on heating than on cooling (inverse hysteresis). The 
conductivity of Rb3D(SeO,)2 in the low temperature phase is 
higher on the first heating than on cooling from the superionic 
phase. In the first heating/cooling cycle of Rb3D(SeO4)2 the 
inverse thermal hysteresis of o, i.e. the temperature difference 
between the parts of o(1/T) curves at the highest slope is more 
than 10K. The hysteresis of Ts, i.e. the temperature difference 
between the low temperature limits of the superionic phase (on 
heating and cooling) is smaller and amounts to ~5 K. On the 
other hand, no hysteresis has been observed in the conductivity 
behavior of Rb3H(SeO4)2. Figs. 3 and 4 show that there is prac- 
tically no difference between the conductivity curves obtained 
on heating and cooling. 

As shown in Fig. 4, the transition to the monoclinic phase on 
cooling is completed at ~450 K in Rb3H(SeO4)2, whereas in the 
deuterated crystal it is completed at ~430 K. The temperatures 
are close to the respective DSC peaks (Figs. 1 and 2). The acti- 
vation energies of ionic conductivity in the superionic phase are 
the same in the limit of the error: Ee =0.49 eV in Rb3H(SeO,)2 
and 0.50 eV in Rb3D(SeO4)2 


3.3. Room temperature Raman spectra and X-ray study 


Fig. 5 shows the NIR Raman spectra of Rb3H(SeO,)2 
and Rb3D(SeO4)2 crystals recorded at room temperature. The 
spectra are very similar but it should be observed that the 
stretching vibrations (DSeO.) appear at higher wavenumber 
(788.9 cm7!) than the w(HSeO,) ones (778.6 cm~!). It is obvi- 
ous, that substitution of a higher mass without any change in 
the H-bond parameters has to result in a decrease in the donor 
frequency. The observed increase in the wavenumber seems to 
be a result of a competition between the higher mass substitution 
and the geometric isotope effect. To learn about the deuteration- 
induced changes in the H(D)-bond parameters, X-ray diffraction 
data of both crystals were collected at room temperature. Table | 
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Fig. 5. Room temperature Raman spectra of Rb3H(SeO4)2 and Rb3D(SeO4)2 
single crystals in the range of skeletal vibrations modes. 


summarizes the basic crystallographic data derived at 295 K and 
Table 2 shows the parameters of the H/D-bond in both crys- 
tals. Apart from the difference in the O-H. - -O and O-D. - -O 
angles, the distance between donor O(2) atom and deuteron in 
Rb3D(SeO4)2 is considerably greater than the donor O(2)—H 
separation in Rb3H(SeOq)2. The increase in w(DSeO4) indi- 
cates therefore, that the effect of deuterium-induced changes in 
the geometry of the dimers prevails over that resulting from the 
mass increase. 


3.4. Temperature evolution of the ferroelastic domain 
structure 


The superionic transition in Rb3H(SeO4)2 and Rb3D(SeO4)2 
is accompanied with a change in the crystal structure from 
the monoclinic 2/m to the rhombohedral R3m. According to 
the symmetry considerations [16,17], such a transformation 
results in a possible appearance of three orientation states in the 
2/m phase. The states correspond to three types of ferroelastic 
domains that can be separated by up to six kinds of admissible 
domain walls. The crystal structure is identical in each of the 
states, but the states differ in the spontaneous strain. The mechan- 
ical stress can transform one orientation state into another, but 
different states can coexist in a crystal with no applied stress 
since they are energetically equivalent. The crystals studied 
transform at Tc into the so-called prototype structure (R3m), 
where the spontaneous strain (and also ferroelastic domain struc- 
ture) vanishes. 

The temperature variation of the ferroelastic domain pat- 
tern shown in Fig. 6, was observed on the (001) surface of 
the virgin crystals. The sample of Rb3D(SeO4)2 was divided 
into two variants at room temperature (one wall separating 
two ferroelastic domains) and Rb3H(SeO,)2 one had a multi- 
domain structure with large single-domain areas. The domain 
pattern did not change on heating up to ~435 K in Rb3D(SeO4)2 
(~22K below Tc) and to ~459K in Rb3H(SeO4)2 (~7K 
below Tc). At these temperatures the thin domains of other 


Table 1 

Crystal data and structure refinement for Rb3D(SeO4)2 and Rb3 H(SeO4)2 
Identification code Rb3D(SeO4)2 
Empirical formula D2 O16 Rb6 Se4 
Formula weight 1088.69 

Temperature (K) 293(2) 

Wavelength (A) 0.71073 


Crystal system, space group 
Unit cell dimensions 


Monoclinic, C2/c 


a= 15.403(3) A, a =90°; b=6.0850(10) A, 
B = 102.95(3)°: c= 10.461(2) A, y=90° 


Rb3H(SeO4)2 

H2 O16 Rb6 Se4 

1086.68 

293(2) 

0.71073 

Monoclinic, C2/c 

a= 15.459(3) A, a =90°; b=6.0910(10) A, 
B = 103.02(3)°; c= 10.482(2) A, y=90° 
961.6(3) 

2, 3.753 


Volume (A?) 955.5(3) 
Z, calculated density (Mg m7?) 2, 3.784 
Absorption coefficient (mm!) 22.933 


Crystal size 
Theta range for data collection (°) 


Reflections collected/unique 
Completeness to theta = 27.47 


Refinement method 
Data/restraints/parameters 
Goodness-of-fit on F? 


Final R indices [J >2o0()] 
R indices (all data) 
Extinction coefficient 


0.26 mm x 0.22 mm x 0.19 mm 
4.27-27.47 


2510/591 [R(int) = 0.0673] 
54.1% 


Full-matrix least-squares on F? 
591/1/65 
1.031 


R1=0.0440, wR2=0.1100 
R1=0.0489, wR2=0.1136 
0.00385(9) 


22.787 
0.31 mm x 0.26 mm x 0.21 mm 
3.61-33.14 


5498/1599 [R(int) =0.1181] 
86.6% 


Full-matrix least-squares on F? 
1599/1/65 
1.167 


R1=0.0541, wR2=0.1199 
R1=0.0711, wR2=0.1255 
0.00235(7) 


Table 2 

Details of hydrogen bonds (A, °) at 295 K (D—donor, A—acceptor, D(1)—deuterium) 

Compound D-H d(D-H) d(H: --A) ZDHA d(D---A) A 

Rb3D(SeO4)2 O(2)-D(1) 0.949 1.621 167.72 2.556 O(2) [—x +1, —y, =z] 
Rb3H(SeO4)2 O(2)-H(1) 0.912 1.676 159.01 2.548 O(2) [—x+ 1, —y, =z] 
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Fig. 6. Temperature evolution of ferroelastic domain structure in a) Rb3H(SeOq)2 and b) Rb3D(SeO4)2. 


orientation states started to appear and their number increased 
on further heating. The changes in the microscope image, 
which we interpreted as the appearance of nuclei of the 
high temperature phase, started at ~450K in Rb3D(SeQgq)2 
(~7K below Tc) and ~464K in Rb3H(SeQO4)2 (~1 K below 
Tc). At their respective Tc, the crystals transformed to the 
paraelastic (superionic) phase and the ferroelastic domains 
disappeared. 


4. Discussion 


The data published so far show slight differences in the 
temperature behavior of Rb3H(SeO4)2 crystal. For instance 
the difference in Ts value, which was reported to range from 
449 K to 456 K [3,7-9] can be related, in our opinion, to dif- 
ferent concentrations of lattice defects in the crystals grown 
in slightly different conditions. Also the number, shape and 
separation of peaks of DSC anomaly in the superionic tran- 
sition region [9,18], as well as the details of the physical 
behavior in the so-called “intermediate phase” (temperature 
evolution of ferroelastic domain structure, symmetry change) 
[9,19] may depend on the conditions in which a given crystal 
was grown. 

Irrespective of the above differences, the proton conductivity 
value and the anisotropy of the activation energy in the supe- 
rionic phase are very similar for various Rb H(SeO4)2 crystals 
[15,20]. In the present work we obtained the activation energy 
along the c-axis as EH = 0.49 eV for Rb3H(SeO4)2 above Ts, 
which agrees with earlier results. 

For Rb3D(SeOq)2, the activation energy was ED = 0.50eV 
in the high temperature phase, which contrasts with E, =0.28 eV 


reported in [10]. We consider EP = 0.50eV to be a reason- 
able value, since the same molecular mechanism is responsible 
for the ionic conductivity in Rb3H(SeO4)2 and Rb3D(SeOq)2. 
Furthermore, both crystals undergo the same change of symme- 
try at Ts and show similar properties (ferroelastic, superionic). 
In the superionic phase of hydrogen sulfates and selenates 
the high dynamics of XO, tetrahedra results in fluctuations 
of O(2) oxygen atoms (donor and acceptor) separation. The 
fluctuations are large enough for proton (deuteron) transfer 
from donor to acceptor in its vibrational ground state, so the 
potential barrier for deuteron transfer is expected to be a lit- 
tle higher (by the energy difference of H* and D* ground 
states) [4]. A possible reason for the discrepancy may be 
the one-frequency method of the conductivity measurement in 
[10]. 

As mentioned, the huge H/D isotope effect in Rb3 H(SeO4)2 
was found in the low temperature range (deuteration-induced 
antiferroelectric phase transition) [12]. On the other hand, we 
observed a rather small difference in the Ts value ae — TP x 
5K) between protonated and deuterated crystals on the first 
heating. The conductivity ratio between Rb3H(SeO4)2 and 
Rb3D(SeO4)2 was approximately 1.6, which points to a classical 
activated diffusion mechanism [13,21]. 

Qualitatively similar high temperature behavior was reported 
for another superionic acid salt, CsHSO, [13,14,22]. The very 
small isotope effect on the superionic proton conductivity in 
Rb3H(SeO4)2 shows that protons (deuterons) follow the dynam- 
ics of their environment without encountering considerable 
barriers along the diffusion path and substitution protons by 
deuterons does not significantly influence the intensive dynam- 
ics of the selenate tetrahedra. 
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Polarization microscope observation of the temperature evo- 
lution of the ferroelastic domain pattern, as well as DSC and o(T) 
measurements indicate, that the temperature range of superionic 
phase transition differs in the two crystals studied. On the other 
hand, the transition in both crystals is accompanied with struc- 
tural reconstruction of the lattice from C2/c to R3m and with 
a great increase in the molecular dynamics. Indeed, the transi- 
tion enthalpy is similar in the two crystals, so one can infer that 
similar energy consuming processes are involved. Our results 
indicate that the processes in Rb3D(SeO4)2 are spread over a 
much wider range of temperatures. 

The temperature evolution of the ferroelastic domain struc- 
ture (Fig. 6) shows that in a virgin crystal new domains start 
to appear in Rb3D(SeO,)2 on the first heating at Tc — T+ 22K 
and in Rb3H(SeO4)2 at Tc -T~ 7 K. Similarly, the high temper- 
ature phase nuclei in Rb3D(SeO4)2 start to grow ~7 K, whereas 
in Rb3H(SeO4)2 only at ~ 1 K below Tc. The energy of a domain 
wall (i.e. the energy of formation of a twin) decreases with 
increasing temperature and disappears at Tc, at which no energy 
is required to form a twin wall. The increasing ease with which 
domain walls can be formed on heating is equivalent to the 
well-known effect of an elastic softening of the crystal struc- 
ture on approaching Tc [23]. The number of walls in a crystal 
depends exponentially on their formation energy. As a result 
a dense domain pattern is observed in both crystals close to Te 
(Fig. 6). The formation of new twins at greater (Tc — T) distance 
in Rb3D(SeOq)2 points to a more elastically “soft” structure of 
the deuterated crystal. 

It is worth noting, that there is a correlation between the 
ferroelastic domain evolution in the two crystals and their 
respective thermal anomalies: the domain pattern starts to 
change on the first heating at a temperature, which is close 
to the position of the respective DSC peak. The temperatures 
at which the nuclei of the high temperature phase start to 
appear (~450 K in Rb3D(SeO4)2 and ~464 K in Rb3 H(SeOq)2) 
are close to those at which the highest slope of o (1/T) is 
observed. Similarly, the temperatures Tc at which the ferroe- 
lastic domain structure vanishes (~457 K and ~465 K) are near 
to the respective Ts values, which set the low temperature lim- 
its of the superionic phase (with Ee ~0.50eV) on the o (1/7) 
dependence. 


5. Conclusions 


DSC, ionic conductivity and ferroelastic domain structure 
studies of the isotope effect on the superionic phase transi- 
tion in Rb3 H(SeO.)2 single crystal have shown that deuteration 
(to a degree of approximately 90%) results in a downward 
shift of the transition temperature in the first heating run, 
but practically does not affect the activation energy value 
(Ec © 0.50 eV) in the superionic phase. The small H/D iso- 
tope effect on the proton conductivity is comparable with 
that reported for CsHSO4 and points to a classical activated 
diffusion mechanism: proton (deuteron) follows the dynam- 
ics of its neighborhood and the dynamics does not change 
significantly on substitution of protons by deuterons. It has 
been found that the temperature variation of the ionic con- 


ductivity in Rb3D(SeOq)2 differs considerably in the first 
heating run from that obtained on cooling and in subsequent 
heating—cooling cycles. Such behavior was not observed in 
the non-deuterated crystal. Furthermore, the ionic conductivity 
study of Rb3D(SeO4)2 has revealed an inverse thermal hystere- 
sis of Ts, whereas the protonated crystal has not exhibited any 
hysteresis. 

The difference in the temperature evolution of the ferroelas- 
tic domain pattern in Rb3D(SeOq)2 and Rb3H(SeO,)2 single 
crystals corresponds to the results of DSC and impedance spec- 
troscopy studies and points to a more diffuse character of 
the superionic phase transition in Rb3D(SeOq4)2. The temper- 
ature range of the transition in Rb3D(SeO4)2 is stretched over 
18-22 K, whereas in Rb3H(SeOx)z2 crystals this range amounts 
only to ~7 K. The room temperature ferroelastic domain pat- 
tern remains static on the first heating up to 7c —-T*¥ 22K 
in Rb3D(SeOq4)2 and Tc — T*~7K in Rb3H(SeO4)2. At these 
temperatures new domain walls start to grow. The greater ease 
of formation of the walls in Rb3D(SeO4)2 indicates that the 
crystal structure is elastically softer compared with that of 
Rb3H(SeO4)2. 

The room temperature NIR Raman spectra of the skele- 
tal modes in Rb3D(SeO4)2 and Rb3H(SeO.)2 crystals point 
to a competition between the conventional isotope effect (i.e. 
decrease in the wave number) and deuterium-induced geometric 
changes in the O3SeO(2)—H(D). - -O(2)SeO3 dimers. 
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